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ABSTRACT 

The  advent  of  the  missile  and  satellite  era  has  brought 
with  it  the  need  for  large  and  highly  accurate  radar  tracking 
systems.   Thus,  the  antenna  tracking  error  generated  by  wind 
disturbances  has  become  a  significant  consideration.   This 
study  presents  a  compensation  technique  which  substantially 
reduces  these  disturbance  errors  without  degradation  of  the 
system  response  to  the  tracking  input. 

The  study  was  proposed  by  Philco  Corporation,  Western 
Development  Laboratories,  Palo  Alto,  California.   The  proposal 
resulted  from  tests  conducted  on  a  production  model  of  a 
satellite  tracking  system  having  a  60  foot,  solid  skin,  high 
inertia  antenna.   A  large  portion  of  this  thesis  was  accom- 
plished at  the  Philco  Laboratories  during  the  Summer  of  1960. 

The  authors  wish  to  express  their  appreciation  to 
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his  assistance  and  counsel  in  the  preparation  of  this  thesis 
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1.   Introduction 

The  purpose  of  this  study  is  to  present  a  compensation 
technique  for  reducing  antenna  tracking  errors  caused  by 
wind-gust  disturbances.   The  antenna  system  upon  which  this 
study  is  based  is  the  position  control  system  for  a  60  foot 
radar  antenna.   This  system  meets  the  desired  response  to  the 
tracking  input.   However,  the  tracking  error  caused  by  wind 
gusting  is  not  acceptable.   The  compensation  utilized  must 
reduce  this  error  without  affecting  the  system  response  to 
the  tracking  input.   This  problem  corresponds  to  having  an 
antenna  system  installed  at  the  tracking  site  which  meets 
the  designed  input  tracking  response.   Data  showed  excessive 
tracking  error  resulting  from  dynamic  wind  pressure  on  the 
solid  skin  of  the  large,  high  inertia  dish.   This  error  must 
be  reduced  by  adding  some  type  of  compensation  that  will  not 
modify  the  designed  tracking  input  response. 

An  error  of  a  few  seconds  of  arc  can  be  excessive  in 
determining  the  precise  position  of  a  high  altitude  orbiting 
satellite.   To  meet  such  stringent  error  regulating  require- 
ments while,  at  the  same  time,  achieving  acceptable  tracking 
input  response  using  conventional  compensation  techniques 
poses  a  formidable  problem.   Compensating  techniques  require 
the  sensing  of  the  effects  of  the  wind  loading  on  the  antenna 
system.   Wind  sensing  might  take  the  form  of  external  devices 
such  as  sensitive  annamometers .   These  might  be  placed  in 
locations  around  the  periphery  of  the  antenna  and  at  distances 
so  as  to  determine  in  advance  the  load  disturbance  which  the 


antenna  would  see.   Although  the  nature  of  a  prediction 
sensing  device  lends  itself  to  the  development  of  the  ideal 
regulation  system,  the  accuracy  of  the  prediction  must  of 
necessity  be  at  least  as  accurate  as  the  tolerance  of  the 
desired  output  response.   To  date,  no  wind  sensing  trans- 
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ducers  have  been  produced  that  measure  accurately  transient 
wind  conditions  over  a  sufficiently  small  increment  of  time. 
Even  if  such  devices  were  feasible  and  available,  their 
sensing  would  not  necessarily  reflect  the  actual  loading 
of  the  antenna  dish  some  distance  away.   A  mountainous  site 
might  be  expected  to  show  large  drafting  and  eddy-currents 
about  the  dish  which  would  preclude  the  accuracy  of  any 
external  sensing  device  to  determine  actual  dish  loading. 

A  second  and  more  reasonable  approach  might  take  the  form 
of  a  sampling  of  the  torque  on  the  antenna  drive  shaft. 
Corrections  would  then  be  made  in  the  control  system  based 
upon  these  samplings.   This  approach  is  complicated  and 
expensive,  and  at  least  part  time  use  of  a  digital  computer 
would  be  required. 

A  simpler  and  more  direct  approach  would  be  to  continuously 
compare  the  output  position  under  external  loading  with  the 
output  obtained  having  no  external  load  with  the  same  tracking 
input.   The  difference  of  these  outputs  would  then  be  a 
measure  of  the  effects  of  transient  wind  loading  on  the 
antenna.   Although  this  approach  in  no  way  predicts  the  loading 
effect,  it  does  permit  the  accurate  measurement  of  the  effect 
which  in  turn  may  be  compensated  at  a  faster  rate  than  that 
possible  by  the  basic  system. 
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The  solution  presented  here  will  be  developed  in  the 
above  manner  using  as  the  basic  approach  the  conditional 
feedback  concept  set  forth  in  the  American  Institute  of 
Electrical  Engineers'  publication,  "Applications  and 
Industry,"  July  1955,  by  G.  Lang  and  J.  M.  Ham.   The  condi- 
tional feedback  system  permits  design  specifications  on 
input-output  response  and  on  disturbance-output  response  to 
be  met  independently. 

To  prevent  repetition  and  to  simplify  the  presentation, 
only  the  azimuth  axis  of  rotation  will  be  considered.   Design 
criteria  for  the  compensation  will  be  developed  based  on 
selected  arbitrary  requirements  for  disturbance-output 
response  improvement.   A  specific  compensation  will  be  designed 
using  the  developed  criteria  and  the  response  will  be 
determined  for  various  disturbance  inputs.   A  design  technique 
is  developed  which  facilitates  the  determination  of  this 
compensation . 

An  investigation  of  the  characteristics  of  a  severe  wind 
gust  disturbance  will  be  utilized  to  develop  a  disturbance 
input  which  will  subject  the  compensated  system  to  the  most 
demanding  conditions  expected. 

An  analog  computer  simulation  will  be  made  of  the  complete 
compensated  system.   Data  will  be  displayed  showing  the 
effects  of  the  selected  compensation  for  various  wind 
disturbance  inputs. 


2.   Development  of  the  basic  approach. 

In  the  conventional  single-loop  feedback  control  system 
the  input-output  response,  as  well  as  the  disturbance-output 
response,  are  dependent  on  the  feedback.   Thus,  in  most  cases 
utilizing  conventional  feedback  control,  a  compromise  of  the 
input  response  must  be  accepted  if  a  satisfactory  disturbance 
response  is  to  be  achieved.   Fig.  2-1  shows  such  a  conventional 
control  system. 


Figure  2-1.   Diagram  of  conventional  control  system. 


In  Fig.  2-1,  G-  is  the  Laplace  transform  of  some  control  system; 
Gp  represents  the  feedback  compensation.   R  is  the  normal 
input  signal  to  the  system;  U  is  the  disturbance  signal  at 
the  output  of  the  system;  C  is  the  final  system  output. 
The  transfer  functions  associated  with  Fig.  2-1  are  readily 
developed  as  follows: 

C  =  G^R  -  G2C)  +  U 


C(l  +  G1G2)  =  GJ[R  +  U 


C  = 


1  ♦  GXG2 


R  ♦ 


1  +  GlG2 


U 


(2-1) 


Thus,  the  feedback  G0  affects  both  input  and  disturbance 


responses.   However,  the  problem  treated  in  this  thesis  is 
one  in  which  the  effects  of  disturbance  must  be  compensated 
without  compromise  of  the  given  input-output  response.   A 
compensation  technique  meeting  these  requirements  has  been 
presented  in  "Conditional  Feedback  Systems"  by  G.  Lang  and 
J.  M.  Ham  [l]  and  will  be  used  as  the  basis  for  the  approach 
of  this  thesis  problem. 

A  basic  configuration  for  a  linear  conditional  feedback 
system  is  shown  in  Fig.  2-2,  where  G..  ,  G_,  R,  Ut  C,  are  defined 
as  in  Fig.  2-1. 
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Figure  2-2.   Basic  conditional  feedback  control  system. 


The  output  C  of  Fig.  2-2  is 

G1(l  +  BG2) 
C  =  1  ♦  G1G2H   R  * 


1  +  GlG2H 


U. 


Let  B  in  equation  (2-2)  be  defined  as 

B  =  G1H. 
Equation  (2-2)  then  becomes 

1 


C  =  Gj^R  + 


1  +  GlG2H 


U. 


(2-2) 


(2-3) 


(2-4) 


Letting  U  and  R  be  alternately  zero, 

f  )u=0  ■  Gi  • 

and 

1        .  (2-6) 


U  )r=0  1    +  G1G2H 

Equation  (2-5)  shows  that,  when  B  is  defined  by  equation 
(2-3),  the  input  response  is  dependent  only  on  G.  .   Thus,  G~H 
can  be  utilized  to  modify  the  disturbance  response  without 
affecting  the  input  response.   The  output  sensing  device 
would  normally  be  described  by  H.   If  it  is  assumed  that 
H  =  1,  then  B  becomes  equal  to  G.  (equation  (2-3)).   The 
purpose  of  the  function  B  is  shown  by  considering  the  system 
operation  in  Fig.  2-2.   The  feedback  signal  (b)  is  equal  to  (a) 
since  H  =  1.   With  B  defined  as  equal  to  G. ,  (d)  is  equal  to 
(a)  with  U  equal  to  zero.   Thus,  the  difference  signal  (e)  is 
equal  to  zero  and  no  compensation  signal  is  developed. 
Clearly,  (d)  is  equivalent  to  the  desired  output  response  and 
thus  B  is  called  a  "reference  model."   The  reference  model 
is  an  undisturbed  representation  of  the  transfer  characteris- 
tics of  the  antenna  drive  system.   This  conditional  feedback 
approach  permits  the  use  of  the  feedback  function  G^  to 
accomplish  the  desired  disturbance  compensation  without 
change  to  the  original  input  response. 


3.   Development  of  basic  system. 

The  basic  system  investigated  in  this  thesis  represents 
the  azimuth  axis  position  control  system  of  an  actual  land- 
based  tracking  radar  having  a  sixty  foot  diameter  antenna. 
The  system  consists  of  a  hydraulic  drive  with  tachometer 
feedback  and  suitable  reduction  gearing,  together  with  cascaded 
electronic  compensation  and  position  feedback  loop. 

Prior  to  any  possible  investigation  of  the  basic  system, 
the  hydraulic  drive  system  transfer  functions  had  to  be  derived. 
This  hydraulic  system  consisted  of  a  servo  control  value  and 
amplifier,  a  D.  C.  motor-driven  constant  displacement  pump  and 
associated  hydraulic  drive  motor,  and  reduction  gearing.   The 
hydraulic  system  was  enclosed  with  a  tachometer  feedback  loop. 

The  block  diagram  of  this  hydraulic  drive  system  is  shown 
in  Fig.  3-1.  _  _  _  —  , 


Compensated  ,„f  l[T" 
irror  signal     l 


Servo  control 

valve 

& 

amplifier 


Hydraulic 
motor 


9. 


h 


Tachometer  (Gf) 


1 

N 


Figure  3-1.   Hydraulic  drive  system. 

The  transfer  function  (G.)  is  defined  as  the  joint  function 
of  the  servo  amplifier,  servo  control  valve,  and  hydraulic 
drive  motor.   Neglecting  the  tachometer  feedback  (Gf),  the 
output  position  angle  to  the  reduction  gearing  is  (9.),  where 

ex  =  GhX  (3-1) 

and  X  is  the  compensated  error  signal  of  the  basic  system. 
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The  block  diagram  of  the  hydraulic  drive  system  with  normal 
antenna  loading  and  external  load  torque  (T  )  due  to  wind  is 
shown  in  Pig.  3-2. 


Figure  3-2.   Hydraulic  system  with  synthesized  loading. 


Letting  Gf  =  0 , 


t=T   ♦  T  =T  +TT, 

1    m     a    m    L  * 


(3-2) 


N 


N 


where  T   is  the  torque  due  to  hydraulic  motor  inertia  (J  ) 
m  M  '  m 

and  is  defined  as 


T  =  J   s  0.  . 

m     m     1 


(3-3) 


The 


load  torque  (T,  ),  as  shown  in  Fig.  3-2,  is 


TT    =   Jrs^0T    +    f T s0T    +    T      =    Kr(9      -    0T  )    . 
L  LL  LL  e  La  L 


©i  =  ©   , 


1   _a 

N 


(3-4) 
(3-5) 


where  N  is  the  gear  reduction  ratio  of  the  hydraulic  system. 
Substituting  equation  (3-5)  into  (3-4)  gives 


TL  =  KL(^1  "  V  =  (JLS   +  fLs)0L  *  Te  ■ 

N 


(3-6) 
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or  (3-7) 

KL  1 

9   =  L 9   _  ± X 

L    N(JLs2  ♦  fLs  +  KL)  1         (JLs2  ♦  fLS  ♦  KL)   e 
Substituting  equation  (3-7)  into  (3-6)  gives 

KL  KL 

TT    =    KT0, = 0,    ♦ 


L        W1   X        N(JLs2    +    fLS    +    V       X         (JLfi2    *    fLS    *    V 

(3-8) 

The  flow  of  oil  out  of  the  control  valve  (Q  )  is  considered 

to  be  divided  into  three  parts,  that  lost  to  system  leakage 

(Q_),  that  lost  because  of  compressibility  of  the  oil  (Q  ), 

and  that  actually  used  to  drive  the  hydraulic  motor  (Q  ). 

Q  =  Qt  +  Q  +  Q  (3-9) 

xv   ^L   xc    xm 

The  flow  out  of  the  control  valve  is  a  function  of  the  input 

error  voltage  (X)  within  some  constant  (K  )  and  is  defined  as 

Q  =  K  X  .  (3-10) 

^v    v 

The  constant  (K  )  is  defined  to  include  the  control  valve 

v 

constant  as  well  as  the  servo  amplifier  gain  constant.   Brown 
and  Campbell  in  "Principals  of  Servomechanisms"  [2]  define 
the  separate  oil  flows  of  the  system  as  follows: 

QL  =  L  P  ,  (3-11) 

where  L  is  the  leakage  coefficient  and  P  is  the  pressure 
differential  of  the  system. 

Qc  =  |  sP  ,  (3-12) 

where  V  is  the  volume  of  oil  under  compression  and  B  is  the 

bulk  modulus   of  the  oil. 


Q   =  dm  sO,  , 
^m         1  ' 


(3-13) 


where  dm  is  the  motor  displacement  and  0-  is  shown  in 

Pig.  3-2.   Dimensions  of  all  constants  are  listed  in  the  table 

of  constants  on  page  12. 

Making  substitutions  of  these  flows  into  equation  (3-9)  gives 


K  X  =  dm  sO.  +  (L  +  ^s)  P. 
v         1        B 


(3-14) 


The  pressure  differential  (P)  of  the  system  is  related  to  the 
torque  output  of  the  motor  (T..)  as  follows: 


T1  =  dm  P  . 


(3-15) 


Substitution  of  equations  (3-15),  (3-3),  (3-4),  into  equation 
(3-2)  gives 


T.  =  dm  P  =  J  s20.  +  -^ 
1  m    1    N2 


-K, 


_(JLs*  +  fLs  +  KL) 


+  1 


0. 


KT 


N(JLs2  +  fLs  +  KL)   e 


(3-16) 


or 


Jn,  S  KT 

P   -*-   Q1   + 


(JLsz  +  fLs) 


dm 


dm  N2   (JLs2  +  fLs  +  KL)    1 


K, 


dm  N  (JLs   +  fLs  +  KL) 


(3-17) 


Substitution  of  equation  (3-14)  gives 
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K  X  =  dm  sO,  +  (L  +  Vs) 
1         B 


J  s2  0t    KT  (Jr  s2  +  fT  s)  0. 

m     1 

dm 


dm  N2(JLs2  +  fLs  +  KL) 


(L  ♦  Js)  KL 


dm  N  (JLs   +  fLs  +  KL) 


T   . 
e 


(3-18) 


If  the  external  torque  (T  )  is  assumed  equal  to  zero,  then, 

y 

solving  for  3—  in  equation  (3-18)  gives 


91    Kv 


0. 


d»  s  *  (1  .  |s)  -^—  *    (L  ♦  |s)  —  2 


(JLs-  ♦  fLs) 


dm 


dm  N*  (JLs  +fLs+KL) 


Since  G.  =  _1   t 


K 


h    s 


'dm  +  (L  -f  Vs)J  s  +  (L  +  Vs)KT 
B   -™-         b   -t 


(JLS  +  V 


dm 


dm  N2  (JLs2  ♦  fLs  ♦  KL) 


(3-19) 
Since  the  system  stiffness  factor  (K,.)  is  relatively  large, 
equation  (3-19)  reduces  to 


BK 


Gh  = 


dm  /  V(Jm  *    JL/N2) 


s   +  (LB  + 


v      LBfT 
)  s  +     L 


dm2  B 


V   N2(J   +  JT)     N2(J  +JT)    V(J   *  JT) 

m     L  m   L        m     L 


N2 


N2 


N2  _ 
(3-20) 
The  friction  coefficient  of  the  antenna  drive  is  considered 
very  small  with  respect  to  the  load  inertia  constant,  thus, 


NT(J   +  J  ) 
m     L 


«  0  , 


N' 
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and  equation  (3-20)  reduces  to 


Gh  = 


N2K  B  dm  /  V(JT  +  N2J  ) 
v  L      m 


s2  +  LBs  +  N2  dm2  B 


V~ 


V(JT  +  N^J  ) 


(3-21) 


The  hydraulic  constants  associated  with  this  system 
are  as  follows: 

N  =  2160 


K   =  0.895 
v 

JL  =  4.2xl06  lb-in-sec2 

J   =  0.09  lb-in-sec2 
m 

B   =  0.27x10   psi 

L   =  60.5xl0~5  in3/sec/psi 

V   =  3.9  in3 


dm  =  0.215  in  /radian 


(referred  to  antenna) 
(referred  to  motor) 


Substituting  these  constants  into  equation  (3-21)  gives 


Gh  = 


1.345x10 


s(sz  +  42s  +  3228) 


(3-22) 


Equation  (3-22)  is  the  transfer  function  of  the 

uncompensated  hydraulic  system  and  was  used  in  the  initial 

development  of  the  complete  basic  system  which  follows.   The 

transfer  functions  of  the  complete  basic  system  chosen  initially 

for  investigation  are  shown  in  Pig.  3-3.   The  values  shown  for 

the  individual  transfer  functions  were  suggested  by  the 

manufacturer  as  representing  the  azimuth  control  system  for 

the  antenna  under  investigation.   The  open-loop  transfer  function 

(F  )  and  the  closed-loop  function  (P  )  for  the  system  of 
o  c 
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Mgh 


24,050(s+1.7)(s+20) 
s(s+10) 


G 


w5?y  11& 

v^        s(s2+42s+ 


450 

3228) 


h 


( .7x.086)s 


2160 


1 

N 


G      =   Cascade    compensation 
c 

G,     =    Hydraulic    drive 

Gf   =    Tachometer    feedback    (tachometer    constant    =    .086) 

N      =   Gear    ratio    (2160) 


Figure  3-3 
Block  diagram  of  initial  basic  system 
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Fig.  3-3  are  as  follows: 

F   =  14.96  x  104(s  +  1.7)(s  +  20)  (3-23) 

°    s2(s  ♦  10)(s2  ♦  42s  +  4038) 

p   s    14.96  +  104(s  +  1.7)(s  +  20)    

c    s5  +  52s4  +  4458s3  +  1.8988xl05s2  +  3.24xl06s  +  5.08xl06 

(3-24) 
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Routh's  cri 

terion 

test : 

1 

4.458xl03 

52 

1.8988x10 

810 

3.14xl06 

-1.172xl04 

5.08xl06 

3.49xl06 

0 

5.08xl06 

0 

4.   Analysis  of  basic  system. 

In  order  to  ascertain  the  stability  of  the  system  of 
Fig.  3-3,  a  Routh's  criterion  test  was  made  of  the  closed- 
loop  transfer  function  (F  )• 

t: 

3.24xl06 

5.08xl06 

0 

0 

0 

0 

The  system  as  given  was  believed  stable,  but  Routh's 
criterion  test  did  not  verify  this  assumption.   Accordingly, 
a  Nyquist  plot  was  drawn  of  the  system  transfer  function 
in  order  to  observe  the  degree  of  instability  and  is  shown 
in  Fig.  4-1.   Within  the  accuracy  of  the  plot,  the  curve 
appeared  to  pass  through  the  Nyquist  point  indicating  marginal 
stability.   For  further  investigation  of  the  system  stability, 
a  Bode  plot,  using  asymptotes  only,  of  F   was  constructed  and 
is  shown  in  Fig.  4-2.   Although  the  plot  crossed  the  zero  db 
line  with  a  negative  six  db  slope,  the  complex  poles  of  F  , 
having  a  natural  frequency  (w  )  of  63.5  rad/sec,  caused  a  hump 
in  the  plot  which,  if  the  magnitude  were  sufficiently  large, 
would  recross  the  zero  db  line  and  result  in  system  instability. 

In  order  to  estimate  the  maximum  amplitude  of  the  hump, 
that  part  of  the  system  transfer  function  (F)  resulting  from 
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the  complex  poles  can  be  written  in  vector  form  as 


F  =  ^ =  MeJN  ,  (4-1) 

(1  -  ^2)+  j2  5^ 
wo       wo 


where:     w   =  break  frequency  of  the  complex  poles 
§  =  damping  coefficient 

M  =  IfI 

N  =  phase  angle  of  F 

The  maximum  absolute  value  of  F  (M    )  is  the  maximum 

max 


w0 


amplitude  of  the  hump.   Letting  w  =  u  and  substituting  into 
equation  (4-1)  gives 

A 

F  =  


o 

(1  -  u  )  +  j2  $  u 


or 


M  =  .  (4-2) 

V(l  -  u2)  +  4  8 2u2 


Differentiating  equation  (4-2)  with  respect  to  u  gives 

dM  =    -2(1    -    u2)(-2u)    +    8  3  2u . 

du        2    [d    -    u2)2   +    4^2u2J   i   [(1    -    u2)2   +    48  2u2] 

(4-3) 

Setting  equation  (4-3)  equal  to  zero  and  solving  for  u  gives 

u  =  (1  -  2  3  2y    .  (4-4) 

Substituting  equation  (4-4)  into  (4-2)  yields 

Chestnut  and  Mayer,  Servomechanisms  and  Regulating 
Systems  Design,  Vol  I,  pp  312,  September  1953. 
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M     = rr-    .  (4-5) 

max    2*  (1  -*  2)> 


Equation  (4-5)  is  an  approximation  of  M    since  it  was 


max 


developed  on  the  basis  that  the  frequency  at  M    was  equal 

r  '  max 

to  the  break  frequency  (w  ).   However,  these  two  frequencies 
differ  slightly. 

The  damping  factor  ( S   )  can  be  computed  from  equation 
(3-23)  as 

2  S  w   =  42  , 

n       ' 

where  w  =  (4038)1  ,   or  *  =  0.33. 
n  ' 


Therefore,  from  equation  (4-5) 


M    =  1.7  =  4.62  db  . 
max 


The  maximum  amplitude  of  the  hump,  4.62  db ,  was  added  to  the 
amplitude  at  the  complex  pole  break  frequency  in  Fig.  4-2  and 
resulted  in  a  system  amplitude  at  this  frequency  of  0.12  db. 

Although  the  above  calculation  is  not  exact,  it  at  least 
shows  the  hump  quite  close  to  the  zero  db  line  and  thus,  the 
system  has  a  very  small  stability  margin  at  best.   A  decrease 
in  system  gain  would  improve  the  stability  by  lowering  the 
hump  below  the  zero  db  line,  but  as  only  5.5  db  can  be  dropped 
before  the  slope  crossing  the  zero  db  line  decreases  to  -12  db, 
this  approach  is  not  satisfactory. 

The  above  marginal  stability  was  a  result  of  the  location 
of  the  complex  poles  of  F  .   These  poles  in  turn  were  placed 
by  the  location  of  the  complex  poles  of  the  hydraulic  drive 
transfer  function  (G.)  and  the  tachometer  feedback  (Gf).   Near 
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cancellation  of  these  latter  poles  was  achieved  by  insertion 
of  a  complex-zero  compensator  (G  )  in  series  with  G^t 
alleviating  the  instability  problem.   The  design  of  the  complex- 
zero  compensator  appears  in  Appendix  I. 

Equation  (3-22)  has  shown  the  complex  poles  of  G,  to  be 
(s   +  42s  +  3228).   The  form  of  the  transfer  function  of  the 
complex-zero  compensator  (G  )  is  shown  by  equation  (1-9)  in 
Appendix  I.   Thus,  the  G   necessary  to  cancel  the  complex  poles 

of  G,  is 
h 

G   =  Wtq  <s2  *  42s  +  3228). 
o    3228 

For  the  purpose  of  this  study,  complete  cancellation  of  the 

complex  poles  of  G,  was  assumed.   In  the  actual  application, 

some  residue  would  remain  since  equation  (1-9)  does  not 

exactly  represent  the  transfer  function  of  the  compensator 

in  Fig.  1-1  of  the  Appendix. 

Placing  G   in  series  with  G,  inside  the  tachometer 
b   o  h 

feedback  loop  results  in  a  new  open-loop  transfer  function, 

n   -  37  (s  +  1.7)(s  +  20)  ,.    ,* 

F   =  = .  (4-6) 

°       s2  (s  +  10) 

Fig.  4-3  shows  the  Bode  diagram  of  the  function  of 
equation  (4-6).   A  comparison  of  this  Bode  with  Fig.  4-2  shows 
that  the  desired  cancellation  has  been  achieved  without  change 
of  the  system  gain  or  location  of  the  remaining  break  points. 
Thus,  the  system  response  remains  the  same  with  the  exception 
of  the  undesirable  effect  of  the  complex  poles. 

At  this  point,  a  limitation  was  placed  on  the  system 
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cascade  compensation  (G  ) .   In  order  that  components  of  less 

stringent  design  might  be  utilized  in  G  ,  the  cascade  gain  of 

2.405x10   was  limited  to  a  value  of  1x10  .   In  order  that  the 

system  gain  of  equation  (4-6)  be  retained,  servo  amplifier 

gain  of  the  hydraulic  system  as  well  as  tachometer  feedback 

gain  were  modified  through  the  following  development. 

The  open-loop  transfer  function  of  the  basic  system  of 

Fig.  4-4  is 

G     G  G. 

Fo  =  W-   (  1  §  G  GhGf  >  •  <4"7> 

o  h  r 

Equating  the  system  gains  of  equation  (4-7)  with  the  desired 
system  gain  of  equation  (4-6)  gives 


K     K  K. 

IT  (  1  *    K  KhKf  >  =  37  •  ^~^ 

o  h  f 

where  K  represents  the  gain  constants  of  the  subscripted 
transfer  functions  of  equation  (4-7)  .   Substituting  K   =  1000 
in  equation  (4-8)  gives 


KoKh         37x2160  =  g0  m 


1  +  K  K.  K.      1000 
o  h  f 

Substituting  Kf  =  0.7x0.086,  where  0.086  is  the  tachometer 
constant,  into  equation  (4-9)  gives 

(1  -  .7x.086)K  K,  =  80 
o  h 

or 

KoKh  =  ~20-95  •  (4-10) 
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Since  all  transfer  function  gains  must  be  positive,  equation 
(4-10)  indicates  that  the  tachometer  gain  must  be  changed  in 
order  to  realize  the  desired  overall  system  gain.   In  order 
that  K  K,  be  a  positive  number,  equation  (4-10)  demands 

(1  -  80Kf)  >    0  .  (4-11) 

Since  the  tachometer  constant  (.086)  can  not  be  changed,  the 
tachometer  gain  must  be  changed  so  that  its  value  is  less 

that  0.1455  in  order  that  K  K,  be  positive.   To  ensure  this 

oh    F 

positive  value,  a  tachometer  gain  of  0.135  was  selected. 


and 


or 


Kf  =  .135x.086 


(1  -  .135x.086x80)K  K.  =  80  ,  (4-12) 

oh' 


K  K,  =  1150 
o  h 


where  K   =  -,^.-Q  ;   therefore 
o    3228  ' 


Kh  =  1150x3228  =  3.72xl06  .  (4-13) 

Equation  (3-21)  defined  K,  as 

N2  B  dm  K 
Kh  =  V(JT  ;  N2*  )  "  1.528x104  kv  .  (4-14) 

L    1N  °m 

Substitution  of  equation  (4-13)  into  (4-14)  gives 

K   =  245 
v 

where  K   was  previously  defined  as  including  the  control 
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valve  constant  as  well  as  the  variable  servo  amplifier  gain 

The  transfer  functions  of  the  basic  system  reflecting 
the  above  changes  are 

3.72xl06 


h    s(s2  ♦  42s  +  3228) 


G   =  Joe  (s2  +  42s  *  3228) 

o     52.ZO 


Gf  =  (  .135x.086)s 

1000  (s  +  1.7)(s  ♦  20) 
c        s(s  +  10) 


The  new  open-loop  transfer  function  (F  )  was  shown  in 
equation  (4-6).  The  new  closed-loop  transfer  function  was 
derived: 

F   =  37  (s  +  1.7)(s  +  20) (  (4_15) 

C    (s  +  1.74)(sZ  +  45.26s  +  724.4) 


This  is  the  basic  system  transfer  function  referred  to  here- 
after as  G.  and  is  shown  in  Fig.  4-4. 

As  was  shown  in  the  development  of  the  basic  approach, 
the  conditional  feedback  compensation  effects  only  the 
disturbance  response  of  the  basic  system.   Fig.  4-5  shows 
the  response  of  the  basic  system  to  a  unit  step  disturbance 
input  (U) ,  where , 

c  =  r^u  . 

o 
The  response  curve  of  Fig.  4-5  was  plotted  using  the  transient 
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response  equation  determined  by  the  graphical  residue 
technique.   The  details  of  this  procedure  are  included  in 
Appendix  II.   The  maximum  overshoot  in  Fig.  4-5  is  10.3% 
and  occurs  at  t  =  0.10  seconds.   Pall-time,  defined  as  the 
time  interval  between  ten  percent  and  90  percent  decay  during 
the  initial  cycle  of  response,  is  0.037  seconds.   The  condi- 
tional feedback  compensation  to  be  attempted  will  try  to  reduce 
significantly  this  fall  time. 

Pig.  4-6  shows  the  frequency  response  of  the  basic 
system  to  a  disturbance  input  of  from  0.5  cps.  to  40  cps. 
It  is  noted  that  as  the  frequency  of  the  disturbance  input 
signal  is  increased,  the  regulating  effect  of  the  basic 
system  deteriorates.   In  order  to  compare  this  frequency 
response  with  that  resulting  after  conditional  feedback 
compensation,  the  system  bandwidth  is  defined  as  that 
frequency  above  which  the  input  disturbance  is  attinuated 
by  less  than  10.62  db ,  or  where  M  is  equal  to  (1  -  .707). 
In  Fig.  4-6  this  value  of  amplitude  is  reached  at  a  frequency 
of  2.10  cycles  per  second.   The  corresponding  phase  angle 
is  (-)78  degrees. 
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5.   Development  of  compensation  technique. 

The  following  is  a  derivation  of  the  conditional 
compensated  disturbance-output  transfer  function 
together  with  a  presentation  of  a  compensation  technique. 
Fig.  5-1  is  a  block  diagram  showing  the  inputs  and  transfer 
blocks  of  the  proposed  conditional  feedback  compensated 

system,  where 

F 

G.  =  t — ; — 7T-  ,  the  closed  loop  transfer  function  of  the 

1        1   +    F  V.  4. 

o    basic  system 
F   =  the  open  loop  transfer  function  of  the  basic  system 
Gp  =  compensator  to  be  determined 
R   =  tracking  input  information 
U   =  load  disturbance  information 
C   =  system  output 


u 

R*                            v+A           a./Os      b    v 

F 
o 

+ 

y. 

G2 

i 

' 

d 

Gi     - 

-^ 

r 

Figure  5-1.   Block  diagram  of  proposed  conditional 
feedback  system. 
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The  transfer  function  of  the  output  (C)  with  respect 
to  R  and  U  is  developed  from  Fig.  5-1.   The  relations  below 
follow  from  Fig.  5-1. 


C=U+c=U+(a-  C)F 

o 


c  =  bF   =  (a  -  C)F 
o  o 

b  =  a  -  C 


d  =  RGj^  -  C 

e  =  RG2Gl  -  CG2 


a  =  R+e  =  R  +  RG2Gi  "  CGo 

Combining  these  relationships  gives 

C  =  U  +  (R  +  RG^G,  -  CG„  -  OF   .  (5-1) 

2  1      2       o 

Rearrangement  of  equation  (5-1)  gives 


it         t1  +  G1G_)F 

1+F   +  F  G„    1+F   +  F  G   *  *  kd-4) 

o  o  2        o     o  2 

F 

Since  1  +  F   =  -zr-  ,    then  equation  (5-2)  becomes 
o    G,  ^ 


U       C1  "**  GlG2)Fo 

C  =  -5 +  £ z;  °  R 

F  F 

7^  +  F  G0    ^  *    F  G. 
G-     o  2    G1     o  2 


1            (1  +  G1G2)FQ 
C  =  p U  +  p " 2  R 

-2(1  ♦  G^)      _2(1  +  Gig2) 


C  =  (1  ♦  F  )(1  +  G.G9)  U  *  G1R  <5"3> 

o        12 
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Therefore,  the  output  due  to  load  disturbance  (U)  is 


[ 


1  ♦  F 


o  J  L 


1  *  G1G2  J 


U 


(5-4) 


The  output  due  to  the  tracking  input  (R)  is 


GjR  . 


(5-5) 


The  input-output  transfer  function  G..  is  the  same  as  that 
for  the  uncompensated  basic  system.   Thus,  this  input-output 
response  remains  unchanged  with  the  addition  of  the  compensation. 

Root-locus  techniques  are  to  be  used  to  select  a  G?  which 
will  give  a  desired  disturbance-output  response.   The  disturbance- 
output  transfer  function  may  be  written  in  the  following  form 
to  facilitate  the  application  of  the  root-locus  technique. 


D  Dn 
o  2 


1  +  F 


1  +  G1G2J 


1  +  F   +  F  G0 
o    o  2 


0  2     o  2    o  2 


(5-6) 


where 


F   = 
o 


N 

_( 

D 


G,  = 


o 


1  +  F 


N 
o 

In  +  d  ) 

o     o 


G2  =  d! 
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The  denominator  of  the  above  function  requires  two  root  loci 
in  order  to  determine  the  final  closed-loop  roots  of  the 
disturbance-output  transfer  function.   Since  it  is  desired 
to  investigate  the  effects  of  several  values  of  the  compensator 
G-,  it  seemed  advantageous  if  a  function,  closely  representing 
the  disturbance-output  function,  could  be  utilized  which 
would  require  only  one  root-locus  plot  to  place  the  roots  of 
the  closed-loop  function.   An  examination  of  the  function 

shows  it  to  be  such  a  function.   This  function  can 


1  ♦  GlG2 


be  rewritten  as 


1       P  D1D0 

♦  Fo =  LJ (5_7) 


1  +  G,G0    1  +  F   +  F  G-    D.O,  ♦  N.N,, 
12         o     o2     12     12 


Nl 
where   G-  =  =r—  * 
i    u1 

The  roots  of  the  function  of  equation  (5-7),  hereafter  referred 
to  as  the  "representing  function,"  are  the  same  as  that  for 
the  desired  disturbance-output  function  of  equation  (5-6). 
Furthermore,  this  function  is  of  a  form  which  requires  the 
use  of  only  one  root  locus  in  order  to  place  the  roots  of  the 
desired  closed-loop  function.   Thus,  if  a  direct  correlation 
between  the  representing  function  and  the  desired  function  can 
be  shown,  the  problem  of  selecting  Gp  through  simple  root- 
locus  techniques  is  made  relatively  easy.   The  development 
of  this  correlation  follows. 
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The  representing  function  of  equation  (5-7)  can  also  be 
written  as 

1  o  2    o  2 


12     o2     o2     o2 


For  the  basic  type  2  system  given,  D   takes  the  form 


s2(r  lS  +  1)(  r  2s  +  1) 


and  represents  the  denominator  of  the  open-loop  transfer 
function  of  the  basic  system. 

Assuming  the  compensator  (G~)  takes  the  form: 

sm(r  s  +  1)(^2S  +  1} 

,  m  ^  O, 

(r  3s  +  i)(r  4s  +  D 

then  the  numerator  of  the  representing  function,  (D  D_  + 

o    2 

N   D~)  ,    is    some   polynomial    in    s    of    the   form 
o    2  ' 

(Asn   +    Bs11-1    +....    +    D)    . 

The  denominator  of  the  representing  function  is  also  of  the 

same  general  form,  (A'sp  +  B'sp~   ♦....+  D').   Since 

no  factored  s  term  exists  in  either  the  numerator  or  denominator, 

the  response  to  a  step  input  of  this  representing  function  will 

include  a  residue  at  the  "input"  pole  inserted  at  the  origin 

and,  therefore,  will  have  a  steady-state  error  constant. 

The  desired  disturbance-output  function  has  been  shown 

in  equation  (5-6)  to  be: 

DoD2 

D  D0  ♦  nV  +  N  N0   »  Which  is  of  the  form: 
o  2     o  2     o  2 
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sz(r  1s  ♦  i)(r  2s  +  D 


.  (r  3s  +  i)(r  4s  +  1) 


A,SP  +  b,p-i  + 


.  .  .  *  D 


Since  this  form  does  not  produce  a  residue  at  the  origin 

of  the  s  plane  due  to  a  step  input,  the  time  response  has 

no  error  constant. 

The  time  response  of  the  representing  function  to  a 

step  input  is  of  the  form: 

-  oc  t 
C,  +  C.e    1   sin(w  t  +  V  ) ,  to  a  second  order 
12  c 

approximation.   As  time  increases,  the  value  of  the  output 
of  the  representing  function  will  tend  to  be  larger  by  an 
amount  C,  tha>n  that  expected  for  the  desired  function  in  which 
C.  is  equal  to  zero. 

If  there  exists  no  phase  difference  (  "2T  )  between  the 
representing  function  and  the  desired  function,  then  to  a 
second  order  approximation,  the  correlation  would  simply  be: 


1  ♦  GlG2  Au(t)  "  Cl  " 


1  +  F 


1  +  GlG2 


A  u(t) 


This,  of  course,  is  not  the  case.   Some  phase  difference 
between  the  two  functions  must  exist  since  they  are  not 
identical.   It  is  this  phase  difference  (  <f  )  between  the 
functions  which  in  addition  to  the  constant  error  term, 
decides  the  similarity  between  the  two  functions.   Some 
correlation  with  respect  to  this  phase  difference  is  necessary 
An  inspection  of  the  desired  transfer  function, 


( 


■)(• 


)  = 


D  D_ 
o  2 


1  +  F  'xl  +  G.G  '    D  D0  +  N  D0  +  N  N0 
o        12      o2     o2     o2 
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shows  that  the  zeros  of  this  function  are  made  up  of  the 

poles  of  F   and  G9 .   An  inspection  of  the  representing  function, 

.  D0(N   ♦  D  ) 

•         x        2   o     o 


12      o  2     o  2     o  2 


shows  this  function  to  have  zeros  that  are  made  up  of  the 
poles  of  G1  and  Gp  .   Both  functions  contain  zeros  which  are 
the  poles  of  G~.   In  addition,  both  functions  have  the  same 
roots.   The  only  difference  between  the  time  response  of  these 
functions  must  be  due  to  the  difference  between  the  placement 

of  the  poles  of  F   and  the  poles  of  G, ,  (N  +  D  ),  the 

^  o  loo 

remaining  zeros  of  the  two  functions.   To  a  second  order 
approximation,  the  phase  difference  (T)    can  be  determined  by 
noting  the  placement  of  the  poles  of  F   and  G1  with  respect 
to  the  complex  pole  of  either  the  representing  or  desired 
function,  since  they  are  identical.   The  phase  difference  (  T  ) , 
determined  using  graphical  residue  methods,  is  as  follows: 

T  =  Sum  of  the  angles  from  the  poles  of  F   to 
the  complex  pole  of  the  representing 
function  minus  the  sum  of  the  angles  from 
the  poles  of  G1  to  the  complex  pole  of 
the  representing  function. 

(5-8) 

Since  the  response  of  the  representing  function  to  a 

unit  step  input,  to  a  second  order  approximation,  is  of 

the  form 

Cl  *  C2e~ait  sin(wct  + Y  )  '  (5-9) 

then  the  response  of  the  desired  function  to  a  unit  step 
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input  must  be, 


C.,e   1  sin(w  t  +  v  +  r  )  . 
3  c 


(5-10) 


The  development  thus  far  permits  the  determination  of 
the  correlating  parameters,  phase  difference  (  ^  )  and  steady- 
state  error  constant  (C.).   Techniques  for  utilizing  these 
parameters  to  predict  the  desired  response  from  the  representing 
response  are  presented  in  the  following  development. 

A  comparison  of  equations  (5-9)  and  (5-10)  show  that 
they  differ  by  a  constant,  a  proportionality  factor,  and  a 
phase  difference  (  ^  ).   The  value  of  C.  can  be  determined  to  a 
second  order  approximation  and  is  shown  below. 


\ 


* 


□  Poles  of  R.F 
Q  Zeros  of  R.F 
X  Poles  of  G1 


® 


□- 


pole  at 
origin 
inserted 
for  step 
input 


Figure  5-2.   Representing  function  pole  and  zero  placement 

The  representing  function,  referred  to  hereafter  in  this 
section  as  R.F.,  has  pole  and  zero  locations  as  shown  in 
Fig.  5-2.   Applying  graphical  residue  techniques  to  Fig.  5-2 
for  a  unit  step  input  gives 


36 


C,  =  Distance  to  zeros  of  R.P.  from  the  origin 

divided  by  the  distance  to  the  complex  poles 
of  R.F.  from  the  origin. 

Since  G~  is  chosen  such  that  R.F.  approximates  a  second  order 

system,  the  zeros  of  R.F.  are  equal  to  the  poles  of  G-,  thus 

2 

(5-11) 


C   = 
Ul 


w 

W  (R.F.) 
_  n 


where  W  (G1  )  is  the  distance  from  the  complex  poles  of  G-  to 
the  origin  and  W  (R.F.)  is  the  distance  from  the  complex  poles 
of  R.F.  to  the  origin.   Equation  (5-11)  shows,  for  example,  that 
if  the  compensated  response  is  to  be  three  times  as  fast  as 

that  of  G, ,  then  W  (R.F.)  should  be  approximately  three  times 

In  ' 

that  of  W  (G1 ) .   In  this  example, 

ci  "  (3)2  -  -11 

Thus,  for  a  significant  improvement  in  response  time,  C-  will 
be  small.   If  the  phase  difference  (  2T  )  is  also  small,  Cp  is 
approximately  equal  to  C-,  since  both  equations  (5-9)  and 
(5-10)  must  equal  unity  at  zero-plus  time.   With  both  C1  and  V 
small  the  response  to  a  unit  step  input  of  the  desired  function 
will  be  similar  to  that  of  the  R.F.   For  small  C.  and  T  the 
overshoot  of  the  desired  function  may  be  predicted.   The 
overshoot  of  the  R.F.  to  a  second  order  approximation  can  be 
determined  from  the  root  locus  using  simple  second  order 
techniques.   The  value  of  C.  is  then  added  to  this  estimate 
of  R.F.  overshoot  and  is  used  as  the  prediction  of  desired 
function  overshoot. 
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Equation  (5-10)  shows  the  phase  angle  of  the  desired 
function  to  be  equal  to  the  sum  of  V    and  Y  .   The  value  of  y 
together  with  its  sign  was  computed  using  equation  (5-8). 
If  a  prediction  of  the  desired  response  based  on  the  R.  F. 
response  is  to  be  made,  some  knowledge  of  y  is  necessary 
in  order  that  the  effect  of  v    ,  the  phase  difference,  may  be 
known.   For  example,  let  y"  be  a  second  quadrant  angle  and  r   be 
positive.   During  the  initial  period  of  response,  the  sine 
term  of  equation  (5-10)  will  take  on  smaller  values  than 
will  the  sine  term  of  equation  (5-9).   This  is  to  say  that 
for  small  C. ,  the  response  of  the  actual  desired  function 
will  be  somewhat  faster  than  that  of  the  R.  F.   Similar 
examples  can  be  given  for  negative  y  or  for  y  in  different 
quadrants.   Suffice  it  to  say  that  a  knowledge  of  the  quadrant 
location  of  v  is  necessary.   The  shaded  area  of  Fig.  5-2  defines 
an  area  where  the  location  of  the  complex  poles  of  R.  F. 
produce  a  y  always  in  the  second  quadrant.   This  area  is 
defined  in  terms  of  known  parameters.   It  is  that  area  which 
lies  to  the  left  of  the  complex  pole  of  G..  ,  where  any  complex 
root  of  R.  F.  will  have  a  damping  factor  ( <f  )  greater  than  or 
equal  to  that  of  G. ,  and  will  have  a  damped  frequency  of 
oscillation  (w  )  greater  than  that  of  the  complex  poles  of 
G-.   This  area,  in  general,  represents  the  possible  locations 
of  the  complex  pole  of  R.  F.  which  will  produce  a  faster 
response  than  G.,  but  will  not  effect  an  increase  in  over- 
shoot since  the  damping  factor  is  at  least  as  great  as  that 
for  G. .   Thus,  the  most  probable  location  for  placement  of  the 
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poles  of  R.F.  lies  in  an  area  in  which  v  is  known  to  be  a 

second  quadrant  angle.   For  complex  poles  of  R.F.  in  this  area, 
the  determination  of  a  positive  *  will  predict  a  faster 
response  and  a  negative  *  will  predict  a  slower  response 
of  the  desired  function  as  compared  to  the  representing 
function  (R.F.).   As  the  absolute  value  of  r    increases,  the 
responses  of  the  two  functions  become  more  dissimilar. 

Accuracy  of  the  prediction  of  the  desired  response 
from  the  representing  response  is  gained  with  experience  in 
using  this  technique.   The  steps  utilizing  the  "representing 
function  technique"  for  determination  of  the  conditional 
feedback  compensation  G^  are  summarized  below. 

1.  Specify  the  desired  disturbance  response  of  the 
compensated  system  with  respect  to  the  original 
basic  system  response. 

2.  Utilize  the  root-locus  method  to  determine  a  G?  that 
will  place  the  closed-loop  poles  of  R.F.  in  a  position 
such  that  desired  specifications  are  approximately  met 
and  the  R.F.  response  is  approximately  of  second  order. 

a.  Place  poles  in  that  defined  area  in  which 

y  is  a  second  quadrant  angle. 

b.  Choose  W   of  R.F.  complex  poles  such  that 

Wn(R.F.) 

w  /~ — ^ —  is  equal  to  the  desired  improvement 
n   1 

in  the  response  time. 

c.  Pick  a  damping  factor  to  meet  desired  overshoot 
requirement . 
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3.   Determine  correlation  parameters. 

2 

*j  n 

a.   Compute  C. 


LWn(R.P.) 


b.   Compute  y  using  rule  of  equation  (5-8). 

4.  Make  first  prediction  of  desired  function  overshoot. 

a.  Using  second  order  approximation  technique, 
estimate  overshoot  of  R.F. 

b.  Sum  the  estimated  R.F.  overshoot  and  C,  . 

5.  Predict  the  desired  function  fall-time. 

a.  Estimate  fall-time  of  R.F.  response  as, 

W  (G.) 

fall-time    (R.F.)    =    fall-time    (G.)    x  w    ;R    p    k 

n 

(5-12) 

b.  Predicted  fall-time  of  the  desired  function 
response  is  the  above  estimate  modified  by  r  . 
Positive  y   predicts  less  fall-time  and  negative r 
predicts  greater  fall-time  than  that  of  the  R.F. 
response.   The  degree  of  modification  is  propor- 
tional to  the  absolute  value  of  7T  . 

6.  Make  second  prediction  of  desired  function  overshoot. 

a.   Increase  first  predicted  value  if  T  is  positive 

and  decrease  if  negative.   The  degree  of  modifica- 
tion is  proportional  to  the  absolute  value  of  V  . 

7.  Compare  predicted  response  with  specified  desired 
disturbance  response.   If  satisfactory,  proceed  to 
verify  response  of  the  actual  system  including  the 
Gp  conditional  feedback  compensator. 
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6.   Determination  of  conditional  feedback  compensator. 

The  representing  function,  -z —  _     ,  will  be  used  to 

12 

determine  the  specific  G_  compensation  through  root-locus 
techniques  in  the  manner  outlined  in  Section  5.   Observation 

of  the  response  of  the  basic  system  to  unit  step  disturbance, 
Fig.  6-1,  shows  it  to  have  the  following  characteristics: 

Fall-time:  .0375  sec. 

Overshoot:  10.3% 

Steady-state  error:      0 

Number  of  oscillations:   1 
In  order  to  effect  a  significant  reduction  in  the  track- 
ing error  of  the  system  due  to  transient  disturbances,  the 
transient  characteristics  of  the  representing  function  were 
arbitrarily  chosen  as  follows: 

Fall-time:  Decrease  by  300%  to  500%. 

Overshoot:  No  increase  over  that  of 

basic  system. 

Steady-state  error:      Small. 

Oscillations:  Essentially  the  same  as 

for  basic  system. 

The  procedure  followed  in  using  the  root-locus  method 

was  to  choose  a  G_  which  would  effect  an  increase  of  W  (G, ) , 

2  n   1  ' 

the  natural  frequency  of  the  basic  system,  by  an  amount  of 
300%  to  500%.   At  the  same  time,  the  real  roots  of  the  basic 
system  were  compensated  in  such  a  way  as  to  minimize  the 
residues  of  the  real  roots  of  the  resulting  representing 
function  (R.F.).   Having  done  this,  the  R.F.  approximated  a 
second  order  system  and  its  time  response  was  generally  known 
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The  simple  correlations  developed  in  Section  5  would  then 
predict  the  desired  function  response. 

The  basic  system  closed-loop  transfer  function  (G..  )  ,  was 
given  in  Section  3  as, 

r    37(s  +  1.7)(s  +  20) 

1  (s  +  1.74)(s  +  22.63  +  jl4.53)(s  +  22.63  -  J14.53) 

therefore , 

j       (s+1.74)(s+22.63+jl4.55)(s+22.63-jl4. 55)  D2 
l-(-G1G2  :  (s+1.74)(s+22. 63+ jl4.55)(s+22. 63- jl4. 55)  D2  ♦  * 

+  37(s+1.7)(s+20)  N2   » 

(6-1) 

N2 
where  G0  =  =—  . 

Several  values  of  G~  were  utilized  in  an  attempt  to 
meet  the  time  response  requirements  arbitrarily  chosen  for 
R.F.   The  following  are  representative  of  the  compensations 
attempted. 

— , — £l_  .  .  .  .  ,  a  lag  function, 
(s+a) 

K(s+b)  .  v      -    r    .. 

(s^a)   •  ■  •  •  »  b  >  a,  a  lag  function, 

-r±-^ — r-^-    .  .  .  .  ,  a  >  b,  a  lead  function, 
(s+a)  ' 

Ks ,  a  derivative  function. 

The  G2  compensation  which  gave  the  desired  position  of 
the  roots  of  the  representing  function  (R.F.)  was, 

r     21.06(s+96)(s+200) 

2  "    (s+20)(s+1000)    ' 

a  lag-lead  compensating  function.   Fig.  6-2  shows  the  root- 
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locus  of  the  complex  roots  of  the  R.F.  utilizing  this  G„ . 
The  complex  roots  are  shown  at 

s   =  -80-J52  ;  s2  =  -80+J52. 

The  real  roots  of  R.F.  were  solved  numerically  and  were  found 

to  be : 

s3  =  -1.712  ;  s4  =  -20  ;  s5  =  -1660  . 

The  natural  frequency  of  the  basic  system,  wn(Gi)>  is  26.9 

rad.  per  sec.   The  natural  frequency  for  R.F.,  W  (R.F.),  is 

95.35  rad.  per  sec.   This  is  an  increase  of  350%  over  W  (G- )  . 
^  n   l 

Section  5  has  shown  the  form  of  the  representing 

80 1 
function  to  be        (^  +  C2e      sin(52t  +y) 

for  the  complex  roots  of  R.F.  determined  above.   The  phase 

angle  ( ^  )  lies  in  the  second  quadrant  since  the  complex  poles 

are  located  in  the  area  that  insures  a  phase  angle  in  the 

second  quadrant.   The  steady-state  error  constant  (C.)  was 

approximated  using  equation  (5-11)  to  be  .079.   Since  the 

roots  of  this  function  are  the  same  as  those  of  the  desired 

function,  the  latter  function,  to  a  second  order  approximation, 

takes  the  form,  C3e"80t  sin(52t  +  H>  +  Y  ) . 

In  order  to  apply  equation  (5-8)  for  determining  the 
phase  difference  (  y ) ,  the  poles  of  F  ,  G. ,  and  the  complex 
poles  of  R.F.  were  plotted  in  Fig.  6-3  and  T  was  measured 
as  +13  degrees. 

The  R.F.,  to  a  second  order  approximation,  should 
give  a  fall-time  of  about  350%  faster  than  that  of  the  basic 
system,  or  applying  equation  (5-12), 

Fall-time  (R.F.)  =  26.9  rad/sec   x  0.0375  sec  =  .01  sec. 

95.35  rad/ sec 
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Since  V  is  positive  and  a  relatively  small  angle,  it. is 

expected  that  the  fall-time  of  the  desired  function  response 
will  be  slightly  less  than  that  estimated  for  the  R.F.,  that 
is,  slightly  less  than  .01  sec. 

The  overshoot  of  the  R.F.,  to  a  second  order  approximation, 
was  0.008  .   This  value  added  to  the  CL  of  0.079  gives  a  first 
prediction  of  desired  function  overshoot  equal  to  8.7%  . 
Since  Y  was  found  to  be  positive,  the  overshoot  of  the  desired 
response  should  be  increased  to  some  value  slightly  greater  than 
that  of  the  first  prediction  of  desired  function  overshoot,  that 
is,  slightly  greater  than  8.77o  .   Thus,  the  predicted  response 
of  the  desired  function  falls  within  the  parameters  selected 
for  the  new  compensated  system. 

The  desired  transfer  function  is, 

.              D       D0(D   +  N  ) 
1 _      o      2   o o_ 

(1+F  )(1+G1G0)    (D  +N  )   CD  D  "+  N  D0  +  N  N-) 
o      12       oo      o2     o2     o2 

The  roots,  (D  D0  +  N  D^  +  N  N_) ,  have  been  found  to  be, 
o  2     o  2     o  2  ' 

(s+1. 712) (s+20)(s+ 1660 ) ( S+80+ j52) ( s+80- j52) . 

Therefore,  the  desired  transfer  function  is, 

s2(s+1.74)(s+10)(s+20)(s+1000)(s+22.63+jl4.55) 

( s+ 1. 74) ( s+ 1. 712) ( s+20)(s+22. 63+ jl4. 55 )(s+22. 63- jl4. 55)   X 

(s+22.63-j!4.55) 

(s+80+j52)(s+80-j52)(s+1660)   ' 

which  after  simplification  is, 

s2(s+10)(s+1000) .         . 

(s+1.712)(s+1660)(s+80+j52)(s+80-j52)   '     ^0_2; 
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Using  the  graphical  residue  method  on  equation  (6-2), 
the  output  due  to  a  load  disturbance  applied  as  a  unit  step 
is, 

C  =  .990e"80tsin(52t  +  2.530)  -  .000965e-1 *712t  + 

0.437e-1660t  . 

Fig.  6-1  shows  a  plot  of  this  response.  The  predicted 

response  characteristics  together  with  those  obtained  from 
Fig.  6-1  are  as  follows: 

Predicted  Actual 

Fall-time:              <  .01  sec  .007  sec 

Overshoot:               8.7%+  8.6% 

Steady-state  error:       0  0 

Oscillations:  --  1 

The  fall-time  has  been  increased  by  535%  over  that  of 
the  basic  uncompensated  system  response.   The  overshoot  has 
been  decreased  by  16.5%  and  no  steady-state  error  results. 
This  improvement  in  the  regulation  response  did  not  result  in 
any  change  in  the  desired  original  input-output  response  of 
the  basic  system. 

In  order  to  check  the  accuracy  of  the  correlation 
technique  rules  developed,  the  response  of  the  representing 
function  to  a  unit  step  input  was  determined  using  graphical 
residue  methods,  and  was  found  to  be, 

C  (R  F  )  =  *0488  +  .709e~80tsin(52t  +  2.308)  + 

.000745e~1*712t  +  .426e"1660t  . 
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Fig.  6-1  shows  a  plot  of  this  R.  F.  response.   The  difference 
in  phase  angles  between  the  R.  F.  and  the  desired  function 
is  equal  to  (2.530-2.308)  radians.   Thus,  the  true  value  of  r 
is  equal  to  0.222  radians  or  12.7  degrees.   The  r   obtained 
through  the  correlation  rule,  equation  (5-8),  was  13  degrees. 
The  estimate  of  the  R.  F.  C.  was  .079,  while  the  actual  value 
of  C.  was  computed  to  be  .0488.   The  accuracy  of  the 
correlation  rules  used  in  the  representing  function  technique 
has  been  shown  to  be  more  than  sufficient.   The  time  saved 
in  using  this  technique  for  determining  the  desired  conditional 
feedback  compensation  was  considerable. 
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7.   Selection  of  disturbance  input. 

The  external  load  torque  on  the  antenna  due  to  wind 
gusting  is  a  function  of  time  with  three  independent 
variables,  wind  velocity,  wind  direction,  and  the  tracking 
direction  of  the  antenna.   To  evaluate  the  compensation,  it 
is  desirable  to  investigate  its  effectiveness  in  reducing 
the  greatest  expected  disturbance  error.   Since  the  three 
variables  affecting  the  disturbance  output  are  independent, 
the  maximum  error  will  occur  when  the  errors  due  to  each 
variable  are  maximum. 

To  investigate  the  effectiveness  of  the  compensation, 
a  wind  velocity  profile  having  the  greatest  expected  severity 

was  used.   Prom  consideration  of  data  pertaining  to  wind 

2 
gust  velocities  ,  a  representative  velocity  profile  was 

developed  which  would  produce  the  maximum  expected  disturbance 

error.   This  velocity  profile  is  shown  in  Fig.  7-1.   A  steady 

25  miles  per  hour  wind  was  assumed  to  exist  prior  to  the 

gust.   The  wind  velocity  then  increased  smoothly  to  a  peak 

of  60  miles  per  hour  in  one-half  second  and  remained  at  this 

value  for  one-half  second.   The  maximum  gust  velocity  normally 

decreases  at  a  slower  rate  than  it  increases  during  the  gust 

build-up.   Since  only  that  part  of  the  gust  profile  resulting 

in  the  largest  disturbance  error  is  of  interest,  the  trailing 

edge  of  the  actual  gust  profile  need  not  be  considered. 

W.  R.  Gregg,  Aeronautical  Meteorology,  pp.  269-274,  1930 
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Note:   The  wind  velocity  profile  was  developed  as 
straight  line  segments  approximating  a 
smoothly  increasing  wind  gust  in  order  to 
utilize  a  function  generator  in  simulating 
the  wind  velocity  input  to  the  analog 
computer  simulation  described  in  Appendix  IV 


Figure  7-1.   Representative  wind  velocity  profile. 
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The  true  wind  direction  was  considered  constant  since 
wind  direction  changes  are  most  likely  to  be  negligible 
during  the  short  interval  of  the  gust  being  used.   Since 
true  wind  has  been  considered  constant  during  the  one-half 
second  of  applied  gusting,  the  change  in  relative  wind 
direction  is  a  function  only  of  antenna  rotation.   The 
targets  tracked  by  this  antenna  are  earth-orbiting  satellites 
which  require  relatively  slow  tracking  rates.   Thus,  the 
change  in  relative  wind  due  to  antenna  rotation  is  negli- 
gible during  the  one-half  second  gusting  interval.   During 
the  gust,  the  antenna  was  considered  to  be  in  that  position, 
relative  to  the  true  wind,  which  would  cause  the  greatest 
disturbance  torque  for  any  given  wind  velocity. 

A  plot  of  wind  moment  coefficient  (C   )  about  the 
F  mgy 

antenna  gimbal  center  vs.  angle  of  wind  incidence  has  been 
prepared  for  the  given  antenna  from  wind  tunnel  tests  made 
on  a  l/30th  scale  model.    The  data  is  reproduced  in  Fig.  7-2 
The  use  of  a  "spoiler,"  fixed  to  the  back  side  of  the  antenna 
around  the  periphery  of  the  dish,  was  used  to  reduce  maximum 
torque  loading  on  the  antenna  by  approximately  50  percent. 

The  moment  about  the  antenna  gimbal  center  (M  )  was 
derived  in  Appendix  III  as 

Mg  =  201'9  CmgVw2   <lb-ft>  •  (7-1) 


3 
Guggenheim  Aeronautical  Laboratory,  California 

Institute  of  Technology,  Report  737,  1960. 
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where  C    =  dimensionless  moment  coefficient  about  gimbal 
£   center 

V   =  wind  velocity  (ft/sec). 

w  '     ' 

For  a  maximum  velocity  of  60  MPH  and  a  maximum  C 

7  mg 

of  0.189  (130   angle  of  wind  incidence),  the  maximum  M 

3 
computed  from  equation  (7-1)  is  137.3x10   lb-ft.   This 

moment,  M  ,  is  the  external  torque  (T  )  applied  to  the 

output  azimuth  drive  shaft. 
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8.   Evaluation  of  the  conditional  feedback  compensation. 

The  evaluation  of  the  designed  conditional  feedback 
compensation  is  made  on  the  basis  of  a  comparison  between 
the  uncompensated  disturbance-output  response  and  the 
compensated  system  response. 

A  step  position  disturbance  input  was  used  in  Section 
6  to  determine  the  design  of  the  compensation  since  this 
input  represented  the  most  severe  disturbance,  though  only 
theoretically  possible.   The  compensation  was  designed  to 
provide  for  an  arbitrarily  selected  improvement  in  system 
response.   These  responses  were  shown  in  Fig.  6-1. 

A  comparison  of  system  response  to  a  step  velocity 
disturbance  input  was  made  to  closer  approximate  the  system 
response  to  wind  gusting.   These  responses  are  shown  in  Pig. 
8-1. 

In  order  to  verify  the  degree  of  response  improvement 
indicated  as  a  result  of  the  above  two  disturbance  inputs, 
the  system  response  to  the  wind  gust  profile  of  Fig.  7-1 
was  obtained.   The  electronic  analog  computer  simulation  for 
this  disturbance  input  was  formulated  in  Appendix  III. 
Appendix  IV  presents  the  analog  simulation  used  to  determine 
the  system  response  to  the  wind  gust.   Both  compensated  and 
uncompensated  system  responses  to  the  gust  profile  of  Fig.  7-1 
are  shown  in  Fig.  8-2a. 

In  order  to  note  the  system  response  to  a  wind  gust 
more  severe  than  that  derived  in  Section  7,  and,  in  addition, 
one  of  less  severity,  wind  gust  profiles  of  one-half  and 
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twice  the  time  interval  of  the  profile  of  Pig.  7-1  were 
applied.   These  responses  are  shown  in  Pig.  8-2b  and  Pig.  8-2c. 
The  frequency  response  curves  for  the  compensated  system  are 
shown  in  Fig.  8-3. 

Table  I  shows  the  disturbance  response  characteristics 
of  the  uncompensated  and  compensated  systems  for  the 
disturbance  inputs  applied.   The  most  significant  parameter 
upon  which  to  evaluate  the  compensated  system  is  the  maximum 
error  caused  by  the  wind  gust.   Table  I  shows  this  error,  due 
to  the  wind  gust  of  Pig.  7-1,  to  be  7.2  seconds  of  arc  for 
the  uncompensated  system.   With  the  addition  of  the  conditional 
feedback  compensation,  this  error  was  reduced  to  0.028  seconds 
of  arc.   It  is  emphasized  that  this  substantial  reduction 
in  antenna  position  error  due  to  wind  loading  has  been  effected 
without  change  to  the  original  input-output  response. 
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9.   Conclusions. 

This  study  has  shown  that  there  exists  a  compensation 
technique  which  can  be  successfully  utilized  to  reduce 
antenna  tracking  errors  due  to  wind  disturbances  without 
effecting  a  change  in  tracking  input-output  response. 

Section  2  introduced  the  conditional  feedback  compensa- 
tion method  and  has  shown  that  this  method  permits  design 
specifications  on  disturbance-output  response  to  be  met 
without  change  to  the  tracking  input-output  response.   A 
specific  conditional  feedback  compensation  was  developed  to 
obtain  a  substantial  degree  of  reduction  of  wind  disturbance 
tracking  error.   The  degree  of  disturbance  error  reduction 
is  presented  in  Table  I  of  Section  8.   The  representative 
wind  velocity  profile  developed  in  Section  7  was  applied  as 
the  disturbance  input,  and  the  resulting  error  decreased  from 
7.2  seconds  of  arc  without  compensation  to  .028  seconds  of 
arc  with  conditional  feedback  compensation.   This  represents 
a  257:1  reduction  in  disturbance  error.   In  addition,  this 
substantial  error  reduction  was  achieved  with  a  faster  response 
time.   The  response  improvement  of  the  compensated  system  is 
made  more  evident  by  noting  that  the  other  response  parameters 
listed  in  Table  I  show  substantial  improvement  over  the  uncom- 
pensated system.   A  significant  disturbance  error  reduction 
might  be  expected  when  it  is  noted  that  bandwidth  has  been 
increased  from  2  cps .  to  98  cps. 

In  order  to  determine  the  system  response  due  to  any 
given  conditional  feedback  compensator,  two  root  locus  plots 
would  have  been  required  to  place  the  roots  of  the  closed  loop 

63 


function.   It  was  anticipated  that  several  values  of 
compensation  would  have  to  be  investigated  before  a  satis- 
factory compensator  could  be  obtained.   The  representing 
function  technique  was  originated  for  the  purpose  of  expediting 
the  above  repetitive  process.   This  technique  was  employed 
in  the  determination  of  the  desired  compensator.   The  represen- 
ting function  technique  proved  to  be  accurate,  expeditious, 
and  simple  to  apply.   A  step-by-step  procedure  for  its 
application  is  given  at  the  conclusion  of  Section  5. 
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APPENDIX  I 
COMPLEX  ZERO  COMPENSATOR 

An  electronic  compensator,  Fig.  1-1,  has  been  developed, 

utilizing  an  operational  amplifier,  which  will  essentially 

cancel  the  complex  poles  of  a  cascaded  transfer  function. 
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Figure  1-1.   Complex  zero  compensator  schematic 
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G.  A.  Korn  and  T.  M.  Korn,  Electronic  Analog  Computers, 
pp.  417-420,  McGraw-Hill  Co.,  1956 
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V2  E   "   V2 

lA   =    Rl  —R^~  =    U  *    15 

l5   "  C3V2S 

I;  =  V  ^  *  *;  +  V> 

V2  "  R3  ♦  R4  ♦  R3R4C3s 

14  =  R0  +  R„  +  R^R.C^s  (1-4) 

3     4     3  4  3 

Substituting  (1-2),  (1-3),  and  (1-4)  into  (1-1)  gives 

Cl"  +  Rl  +  r/,  RlR2C2s  *  K3  +  R4  ♦  R3R4C3S  =  °  "^^ 

Letting  R1  =  R?  =  R-  =  R.  =  1  megohm,  then  equation  (1-5) 
becomes 

C  es  ♦  r yo   +  ft " To   =    °     (I~6) 

2x10   +  C2sxl0      2x10°  +  C3sxl0  z 

e  (2xl06C1s  +  C1C2s2xl012  +  1)         _£ 


2xl06  +  C-,sxl012  2xl06  ♦  C,sxl012 

Z  J        (1-7) 

If  C-  is  made  equal  to  C-,,  the  denominators  of  (1-7)  become 

equal,  and  (1-7)  reduces  to 

-  =  -C1Coxl012(s2  +  2  v1a6  s  +  i ™  ) 

e     12  C2xl0o      C1C2xl012  ;  " 

(1-8) 
If  the  capacitors  are  valued  in  micro-farads,  equation 
(1-8)  becomes 

|  =  -C1C0(s2  +  S-  s  +  ^t-tt-    )  . 

e      X  d  U2     U1U2  (1-9) 
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Equation  (1-9)  is  of  the  form  of  a  constant  times  a  complex 
zero,  where  C,  and  C^  are  picked  so  that  the  coefficients 
in  this  equation  equal  those  of  the  complex  poles  it  is  to 
cancel . 

Depending  on  the  values  of  capacitors  used,  it  may  be 
necessary  to  modify  the  circuit  values  slightly  to  suppress 
high  frequency  oscillations  in  the  compensator  amplifier. 
Reduction  of  the  value  of  C~  or  placing  a  relatively  small 
resistance  (about  0.1  megolms)  between  C,  and  ground  has 
been  found  effective.   Such  changes  naturally  reduce  the 
cancellation  effectiveness  of  the  compensator,  but  to  a 
degree  generally  within  acceptable  limits. 
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APPENDIX  II 

GRAPHICAL  RESIDUE  TECHNIQUE  FOR  DETERMINING 
THE  TRANSIENT  RESPONSE  EQUATION  OP  A  SYSTEM 

Consider  the  disturbance-output  transfer  function 
of  the  uncompensated  basic  system: 

C  _1 _     s2(s   +    10) 

U         1    +    F  (s    +    1.74)(s    +    22.63    -    jl4. 55)      X 

o 

X    (s   +    22.63   +    jl4.f!5)     *  CII-D 

Let  U(t)  =  A  u(t),  where  A  =  1.   Therefore  U(s)  =  -  ,  and 
equation  (II-l)  becomes 

r  .  s(s  +10)     

(s  ♦  1.74)(s  +  22.63  -  jl4.55)(s  +  22.63  ♦  J14. 55)  * 

(II-2) 

The  poles  and  zeros  of  equation  (II-2)  are  plotted  in  Fig.  II-l 

The  residue  at  (-1.74)  is  equal  to 

T(  vectors  from  the  zeros  to  the  pole  )  e  P° 

P (  vectors  from  the  remaining  poles  to  the  pole) 

or,        (r1  -  z1)(r2  -  z^)  r-t 

[>!-<«+  Jwe)]  [r,  -    (  a  .  j„c)]  e    •   (II3) 

From   Fig.    II-l    equation    (II-3)    becomes 

(1.74    /l80  °  )(8.26     /0°   )  1.74     /l80°    t 


(25.45  /2  TT   -    Q1)(25.45    /&.) 


-1    74    t 
-0.0222    e      -1--  /i+    l   .  (II-4) 
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The  residue  at  the  complex  pole  (  -  22.63  +  J14.55) 
is  equal  to: 

[(-   oc  +    jwc)    -    z1        [(-    a  +    jwc)    -    z2  (-  oc   +    jw    )t 


c 
e 


R-    a  +    jvvc)    -    r±        Q-   oc  +    jwc)    -    (-   oc  _    jwc)J 

(11-5) 

From  Fig.  II-2,  equation  (II-4)  becomes 

(26.9  /I.  X19.25  /09)  ,  _0  ,_    .,.  „.. 

Z 2.  (-22.63  +  jl4.55)t 

=  e 

(25.45/TT-Q1)(2  x  14.55   /Q-) 

(26.9  ^V  )(19.25  ^/Q^)         e  (-22.63  +  jl4.55)t 

=  x  

(  25.  45/TT-Q1)(  14.  55)  2j 

=  1.395  e  "22'63tsin(14.55t  +  T  ♦  92  -  TT+  c^) 

=  1.395  e  "22,63tsin(14.55t  +  2.319)  .  (II-6) 

Combining  equations  (II-6)  and  (II-4)  yields  the  complete 

transient  response  equation  for  the  system    1     u(t)  : 

1  +  F 
o 

C  =  1.395  e  "22'63tsin(14.55t  +  2.319)  -  0.0222  e  "1-74t 

By  this  same  technique,  the  other  disturbance-output 
transfer  functions  used  in  this  thesis  can  be  converted  to 
their  equivalent  transient  response  equations  for  various 
input  signals.   A  listing  of  these  transient  response 
equations  for  unit  step  and  unit  velocity  inputs  follows. 
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C  1 _  s^(s   +10)  ^__^ 

U   "    1    +    F  (s    +    1.74)(s    +    22.63    -    jl4.55)(s    +    22.63    +    jl4. 55) 

t^-ttt   =    1-395    e    "22-63t    Sin(14.55t    +    2.319)    -    0.0222    e    _1-74t 
lu(  t) 


CJ}\    =    0.0520    e    "22-63t    sin(14.55t    -    0.247)    ♦    0.0127    e    _1-74t 
tu(  t) 


C 

U 


1    ♦    GlG2 


(s  +  1.74)(s  +  22.63  -  j!4.55)(s  +  22.63  +  j!4.55)(s  +  1000) 
(s  +  1.412)(s  +  80  -  j52)(s  +80  +  j52)(s  +1660) 


?(Hv  =    0.0488  +  0.709  e  "80t  sin(52t  +  2.308)  +  .000745  e  ~1'712t 
lu(  t) 


+  0.426  e 


-1660t 


t^ttt  =  0.00172  +  0.0488t  +  0.00738  e  "80t  sin(52t  -  0.265) 
tu(t) 

+  0.000423  e  ~1-712t    _  0.000257  e  -1660t 


C 

u 


1  +  F 
o 


1  +  GlG2 


s  (s  +  10)(s  +  1000) 

(s  +  1.712)(s  +  80  -  j52)(s  +  80  +  j52)(s  +  1660) 


^u(^)  =  °-"°  e  ~80t  sin(52t  +  2.530)  -  0.000965  e  ~1-712t 


+  0.437  e 


-1660t 


tu(t)  =  °*0104  e  ~80t  sin(52t  -  0.0306)  +  0.000564  e  ~1-712t 


-  0.000263  e 


-1660t 
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APPENDIX  III 
DEVELOPMENT  OF  DISTURBANCE  TRANSFER  FUNCTION 

The  purpose  of  this  Appendix  is  to  develop  the  relationship 

between  the  disturbance  input  (U),  a  function  of  the  disturbing 

wind  velocity  (V  ),  and  the  antenna  position  angle  (9r ) .   The 

W  L 

development  of  this  transfer  function  was  undertaken  in  two 

steps.   First,  the  transfer  equation  of  the  disturbance  torque 

(T  )  with  respect  to  V   was  determined.   The  transfer  function 
e  w 

of  T   to  the  antenna  position  angle  (0T )  was  then  determined 
and  the  two  steps  combined  to  give  the  transfer  function  of  V 


w 
to  0L. 

Fig.  7-2  presents  a  plot  of  the  non-dimensional  moment 

coefficient  (C   )  of  the  antenna  as  a  function  of  relative 
mg 

wind  direction.   Using  Fig.  7-2  to  determine  C   ,  the  moment 

b         &  mg ' 

about  the  antenna  gimble  center  (M  )  which  is  also  the  external 


torque  (T  ),  was  found  as  follows: 


T=M=CqSr  (III-l) 

e     g    mg  4  KX±±      * 

O  1  "K 

where,     q  =  ■§- ?  V   ,  the  dynamic  pressure  (t-t2), 

2 
S  =  2830  ft  ,   antenna  dish  frontal  area, 

f   =  air  density  (lb""j?ec2)  , 

ft4 

r  =  30  ft,  antenna  dish  radius. 

The  data  of  Fig.  7-2  was  taken  at  a  temperature  of  20°C . 

and  a  barometric  pressure  of  742  mm  Hg . 

The  air  density  ( ?   )  is  defined  as,    f    -   -%=  /TTT  o\ 

gKl      (.111-2; 

^J.  H.  Dwinell,  Principles  of  Aerodynamics,  pp.  25-26, 
51-52,  McGraw-Hill  Co.,  1949 
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where,     p  =  absolute  pressure  ( — rr) 

ft2 

g  =  32.2  p  ,  acceleration  due  to  gravity 

sec 

R  =  53.3  —  ,  gas  constant 
°R 

T  =  absolute  temperature  (  R)  . 

The  temperature  of  28°C  is  equal  to  an  absolute  temperature 

of  541. 9°R.   A  pressure  of  742  mm  Hg .  is  equal  to  2075  lbs/ft2 

at  28°C.   Substituting  these  constants  into  equation  (III-2) 

gives, 

?    =  .002235  lb-sec2/ft4. 

Solving  for  q  in  terms  of  V   and  substituting  into 

equation  (III-l)  gives, 

T   =  201.9  C    V  2  ft-lbs  .  (III-3) 

e  mg   w 

Equation  (3-18)  of  Section  3  shows  the  relationship 
between  the  error  signal  (X)  and  the  external  load  torque 
(T  )  with  respect  to  the  position  output  angle  (9.)  of  the 
hydraulic  motor.   Fig.  3-2  shows  the  hydraulic  system  configu- 
ration.  Thus,  the  actual  input  to  the  hydraulic  function  (G.) 
f  Fig.  3-2  for  an  applied  external  wind  torque  is  the  sum 
f  the  error  signal  (X)  and  some  function  of  external  torque. 
The  hydraulic  motor  output  angle  (0,)  can  then  be  written  as, 

©!  =  Gh  X  +  G±Gh   Te  (III-4) 

where  G1  is  some  transfer  function  which  carries  external 
torque  into  an  error  position  to  be  applied  to  the  hydraulic 
transfer  function  (G.).   Letting  X  in  equation  (III-4)  be  equal 
to  zero,  then    Q1  =  G^  Tq    .  (III-5) 


o 


o 
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Substituting  X  equal  to  zero  into  equation  (3-18)  gives, 


0. 


(-)  KL  (L+Xs) 


Te    dm  N(JLs2+fLs+KL) 


dm   s+(L+^s)J    s2+(L+^s)KT       (JTs2+fTs) 
B  m  B  L    _       L  L 


dm 


dmN2(JLS2+fLs+KL) 


(III-6) 


Substituting  equation  (III-6)  into  (III-5)  gives, 


GlGh  = 


(-)  KL  (L+|s) 


[dm2  N(JLs2-ffLs+KL)  +  (L+^s)NJins(JLs2*fLs+KL)  +  (L+^s)\(JLs+fL)~] 


N 


(III-7) 


Making  the  same  assumptions  as  in  Section  3,  that  is,  Kr  is 

fL 

very  large  and  —5 =  0,  equation  (III-7)  simplifies  to, 

NZ(J  +JT) 

m   L 

N2 


GlGh  = 


(-)  (L+|s)NB/V(N2Jm+JL) 


s2+BLs+B  dm2N2    ~l 

v    V(N2J  +JT ) 
m   L 


(III-8) 


Dividing  equation  (III-8)  by  G,  as  given  in  equation  (3-21) 


gives, 


Gl  = 


(GxGh)    -(L*|s) 

GT    =  NK   dm" 
h        v 


(III-9) 


Fig.  III-l  shows  the  final  block  diagram  of  the  hydraulic 
system  with  the  added  complex  zero  compensation  (G  )  and 
tachometer  feedback  (G^)  described  in  Section  4.   In  addition, 
the  application  of  T   through  the  synthetic  transfer  function 
G1  is  represented. 
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— KS- 


^> 


9. 


>  1/N 


0, 


Figure  III-l.   Block  diagram  of  final  hydraulic  system 

Assuming  that  0..  is  due  to  the  external  torque  alone, 
(ie.  X'  =  0) ,  then 


9, 


GlGh 


Te    N(l  +  GfGoGh)   • 


(111-10) 


where,     Gr  =  K^s,  the  tachometer  feedback 

K  x  (complex  roots  of  G,  ), 
compensator  added  in  Section  4. 


G   = 
o 


K   x  (complex  roots  of  G,),  the  complex  zero 


Substituting  equation  (III-8)  into  (111-10)  and  simplifying 
gives, 


0L      (s+BL/V)  /  (N  Jm+JL) 


e    s 


~s2+BL/Vs+   B  dm2  N2 


V(N^J  +JT) 
m   L 


!+  KfK0KvBNz  dm 

V(JT+N^J  ) 
L  m 


(III-ll) 
Substituting  the  hydraulic  flow  coefficients  listed  in 
Section  3  and  the  gain  constants  developed  in  Section  4  into 
equation  (III-ll)  gives, 


0 


,-8 


WL  _  1.51x10  °(s+42)     r ad/lb- in 

~ 5 

e   s(s  +42S+3228) 


(111-12) 
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Combining  equations  (111-12)  and  (III-3)  enables  the 
solution  of  the  output  position  angle  of  the  antenna  for 
any  given  wind  velocity.   Thus,  for  a  constant  C   , 

°L   _  201.9  Cmg  x  1.51xl0~8(s+42)  (111-13) 

Vw2      s(s2  +  42s  +  3228) 
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APPENDIX  IV 

ANALOG  SIMULATION  OF  CONDITIONAL 
FEEDBACK  COMPENSATED  SYSTEM 

Fig.  5-2  shows  the  block  diagram  of  the  complete 
compensated  system  under  study.   It  has  been  shown  that  only 
the  response  of  the  system  to  a  disturbance  input  is  of 
interest  since  the  response  to  the  tracking  input  (R)  is 
known  to  remain  unchanged  with  the  addition  of  the  conditional 
feedback  compensation.   Fig.  IV-1  shows  the  block  diagram  of 
the  compensated  system  with  the  tracking  input  equal  to  zero. 


Figure  IV-1.   System  with  conditional  feedback  and  R  =  0. 

The  transfer  blocks  of  Fig.  IV-1  were  defined  as  follows 
in  Section  4: 


o 


_  37(s  +  1.7)(s  +  20) 
s2  (s  +  10) 


r   =  21.06(s  +  96)(s  +  200) 
b2        (s  +  20)(s  +  1000)  ' 

and  where  U  is  the  disturbance  input,  (Or ) ,  the  antenna 

position  due  to  an  applied  wind  velocity.   Equation  (111-13) 

has  shown  the  transfer  function  from  wind  velocity  (V  )  to 

J         w 
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0T  to  be , 

U    o^_    201.9  Cmg  x  1.51  x  10~8(s  +  42) 
Vw2    Vw2      s(s2  +  42s  +  3228) 

The  analog  simulation  of  the  compensated  system  of 
Fig.  IV-1  was  undertaken  in  order  to  determine  the  effects 
of  the  added  conditional  feedback  compensation  on  the  system 
response  to  a  simulated  wind  velocity  disturbance  profile. 

The  following  is  a  list  of  equipment  utilized  in  the 

simulation  of  the  compensated  system: 

3-  Donner  analog  computers,  Mod.  3000,  10  channels 
per  computer. 

1-  Brush  dual-channel  D.C.  amplifier. 

1-  Brush  dual-channel  pen  recorder,  Mod.  BL-202. 

1-  Donner  variable-base  function  generator,  Mod. 
3750. 

1-  Donner  electronic  function  multiplier,  Mod.  3730. 

1-  Donner  null,  voltmeter,  Mod.  5002. 

1-  Hewlett  Packard  low  frequency  function  generator, 
Mod.  202A. 

1-  Moseley  Autograph  x-y  recorder,  Mod.  2-S. 

The  simulation  of  the  disturbance  input  function  was 

undertaken  in  two  parts.   First,  the  disturbance  torque  (T  ) 

was  simulated  using  equation  (III-3), 

T   =  201.9  C    V  2   ft-lbs 
e  mg   w 

The  representative  wind  velocity  profile  of  Fig.  7-1 
was  simulated  as  a  function  of  time  using  a  variable-base 
function  generator.   The  simulated  profile  was  generated  using 
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a  step  velocity  input  as  a  time  base  and  an  eleven  straight- 
line  segment  approximation  of  the  wind  velocity  profile.   The 
constant  initial  velocity  of  25  MPH  was  then  added  to  the 
output  of  the  function  generator.   The  resultant  was  the 
simulated  profile  of  Fig.  7-1.   This  velocity  profile  was 
squared  using  the  function  multiplier.   Since  C    was  considered 

MO  r  mg 

constant  during  the  interval  of  the  gust,  the  squared  velocity 
was  multiplied  by  the  appropriate  constant  to  give  T  . 

Equation  (111-12)  gives  the  transfer  function  between 
T   and  the  output  position  due  to  wind  loading  (9T ) .   The 
analog  simulation  of  this  transfer  function  requires  an 
integrating  output  stage.   This  equation  represents  an  open- 
loop  transfer  function.   Therefore,  there  is  no  feedback  loop 
around  this  output  stage.   Since  the  input  to  this  integrator 
is  non-zero,  this  integrating  stage  would  saturate.   The 
simulation  was  modified  to  include  this  integrating  stage 
within  the  feedback  loop  of  the  servo  system,  resulting 
in  an  output  from  the  disturbance  simulator  of  the  velocity 
(er ) .   Modifying  equation  (111-12)  gives, 

U     ©L    K_  (s  +  42) 


e 
T    T   "  —^ 


e     e    (s*  +  42s  +  3228) 


—  8 
where  K   is  equal  to  1.51x10 
e     M 


Th 


e  velocity  0L  was  added  to  the  velocity  output  of  the 


transfer  function  F   and  integrated  to  obtain  the  system 
position  output  (C)  . 

Fig.  IV-2  shows  the  complete  schematic  of  the  analog 
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computer  simulation  used.   Values  of  components  used  as  well 
as  the  primary  outputs  are  indicated  on  the  figure. 
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